Formation of a Propeller Structure by a Moonlet in a Dense 

Planetary Ring 



Shugo Michikoshi 1 and Eiichiro Kokubo 1 ' 2 
michikoshi@cfca.jp and kokubo@th.nao.ac.jp 

ABSTRACT 

The Cassini spacecraft discovered a propeller-shaped structure in Saturn's A 
ring. This propeller structure is thought to be formed by gravitational scattering 
of ring particles by an unseen embedded moonlet. Self-gravity wakes are prevalent 
in dense rings due to gravitational instability. Strong gravitational wakes affect 
the propeller structure. Here, we derive the condition for formation of a propeller 
structure by a moonlet embedded in a dense ring with gravitational wakes. We 
find that a propeller structure is formed when the wavelength of the gravitational 
wakes is smaller than the Hill radius of the moonlet. We confirm this formation 
condition by performing numerical simulations. This condition is consistent with 
observations of propeller structures in Saturn's A ring. 

Subject headings: Methods: numerical, Planets and satellites: rings 



Introduction 



A moonlet embedded in a planetary ring t ends to open a gap , such as a Keeler or 



Encke gap, due to gravitational scattering (e.g.. iLissauer et al.lll98ll ). Conversely, viscous 
diffusion of ring particles tends to close a gap. If a moonlet is sufficiently large, it will form 
a fully circular gap, whereas a small moonlet will form only a partial gap that consists of 
two azimuthally aligned lobes shaped like a propeller (hereafter, we ref er to this structure as 



a prop eller). Using the viscous fluid model. ISpahn fc Sremcevid (120001 ) and ISremcevic et al. 
(120021 ) predicted the formation of prop ellers by small moonl ets. The Cassini spacecraft 



discovered propellers in Saturn's A ring (ITiscareno et al.l 120061 ) . The moonlets are so small 
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that they cannot be directly detected. Since the radial separation between the two lobes is 
related to the Hill radius of a moonlet, the size of a moonlet with a propeller can be estimated. 
Most of the known propellers are concentrated within narrow bands i n the A ring. Moonlet s 
with propellers have radii ranging from tens of meters to a kilometer (ITiscareno et al.ll2010l ). 
The size distribution o f moonlets has a steeper power-law index than that of ring particles 
( ITiscareno et al.ll2008l ). 



An A^-body simulation is a powerful tool for s tudying ring dyn amic s where collisions 
and g ravitational interactions play important roles. ISeiB et al.l (120051) andlLewis fc Stewart 



(120091 ) performed A^-body simulations of propeller formations. ISeiB et al. considered 
collisions between particles but did not take into account of the self-gra yity of particles. They 
treat ed the effect of self-gravity as increasing the vertical frequency (jWisdom &: Tremaine 
19881 ). They confirmed th e formation of pr o peller s and investigated the sca ling law of the 
propeller size discussed in ISremcevic et al.l (120021 ) . iLewis fc Stewartl (120091 ) included self- 
gravity and the size distribution of particles. They adopted a low optical depth r = 0.1 — 0.2, 
which is smaller than th at of the A ring. Self- gravity causes spontaneous formation of 



gravitational wakes (e.g., ISalol Il995t iDaisaka fc Idalll999l ). They showed that self-gravity 
with a size distribution destroys the propeller structure when the ratio of the moonlet mass 
to the upper limit of the mass distribution is less than 30. 

In this study, we perform real-scale A^-body simulations of a small moonlet embedded 
in dense rings. The optical depth of the A ring is as high as about 0.3-0.5, while that of the 
B ring is larger than unity. Distinct and large gravit ational wakes form in rings with large 
optical depths (e.g., ISalolll995t IDaisaka fc Ida! 119991 ) . Such large gravitational wakes may 
alter the structures around the embedded moonlet. We investigate the condition for propeller 
formation in a dense ring in which gravitational wakes are prevalent. Section 2 summarizes 
the calculation method. In Section 3, we present the simulation results and derive the 
condition for propeller formation. Section 4 consists of a discussion and a summary. 



2. Calculation Method 

We consider a small computational domain of a ring such that L x ,L y <C a, where L x and 
L y are the width and length of the computational domain and a is the distance from Saturn. 
We introduce a local Cartesian coordinate system (x,y,z). The origin revolves around 
Saturn with the Kepler angular velocity Q and is located at the center of the computational 
domain. The x-axis is directed radially outward, the y-axis is parallel to the direction of 
rotation, and the z-axis is normal to the x-y plane. The co mputational domain has periodic 



boundary conditions and is surrounded by eight copies (e.g.. [Wisdom fc Tremaindll988l ; ISalo 
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1995 1) . The equation of motion is linearized, which is referred to as the Hill equation (e.g., 
Petit fe Henonl Il986l ; iNakazawa &: Idalll988l ). The equation of motion is integrated using 
a second-order leapfrog scheme with a variable time step. The sizes of the computational 
domain are L x = 10rn and L y = 80rn, where 7h is the Hill radius of the moonlet given by 
rn = (M/3M s ) 1//3 a, where M and M s are the masses of the moonlet and Saturn, respectively. 

We consider the self-gravity of ring particles and calculate it by directly summing the 
gravitational interactions of all pairs. We calculate not only the gravit ational interactions 
inside the domain but also those from the surrounding copies (jSalolll995l ). The gravitational 
interactions, which are the most computationally expensive part of iV-body si mulations, are 



calcu lated using a programmable, special-purpose computer, GRAPE-DR (IMakino et al. 
2007h . 



The frictionless impact model for hard spheres is adopted (e.g.. lRichardsonlll994l ; lDaisaka fc Ida 



19991 ). In a collision, the normal component of the relative velocity is reduced by a factor e, 
which is the restitution coefficient in the normal direction. The tangential component of the 
relative velocity is conserved. We adopt a restitution coefficient m odel that was deter mined 
by a laboratory experiment. The normal restitution coefficient is ([Bridges et al.lll984j ) 



0.34 min 



where v n is the normal impact velocity. 



lcm/s 



-0.234 



We consider the size distribution of particles. We adopt a power-law model ndR = 
(R/R )~ q dR for -R m i n < R < R max , where R is the radius of particles, R is the radius for 
normalization, q is the power-law index for the distribution, and R m i n and R m ™ are the 
minim um and maximum sizes of particles, respectively. We adopt q = 2.8 ( Zebker et al. 



19851 ). The den sity of particle s is 0.5g/cm 3 . The initial Toomre parameter Q = f2c/3.36G£ 
is set to Q = 2 (lToomrdll964l ). where c is the velocity dispersion of ring particles. 



The moonlet is fixed at the origin of the coordinates. We assume that the moonlet 
radius is 150 m and its bulk density is 0.9g/cm 3 . The semimajor axis and the Hill radius of 
the moonlet are a=117000 km and rn = 228 m, repectively. This semimajor axis corresponds 
to the outer region of the B ring. 
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Results 



3.1. Formation and Non- Formation of Propellers 

We demonstrate formation and non-formation of propellers with low and high surface 
density ring mo dels. The typical m inimum surface density in the B ring is estimated to be 
240-480 g/cm 2 (IRobbins et al .11201 On, whil e the surface density at its outer edge is estimated 
to be 30-70 g/cm 2 (ISpitale fc Porccj UQM) ■ We adopt S = 62 g/cm 2 and 414 g/cm 2 for the 
low and high surface density models, respectively. The size range of particles is in ferred from 
stellar occultation as -R m i n — 30 cm and -R max — 20 m ([French &: Nicholson! |2 00 0l ). However, 
due to the limitations of the available computing resources, we have to adopt the larger 
and smaller values respectively for -R m i n and R m3jX of -R m i n = 2 m and R m -^ = 10 m to fit the 



surface density to the value inferred from the density wave (e.g.. iTiscareno et al.ll2007l ). This 
small size range gives somewhat unrealistic dynamical optical depths of r = 0.18 and 1.2 for 
the low and high surface density models, respectively. However, the surface density controls 
the basic dynamics of dense rings in which gravitational wakes develop. As is shown later, 
the condition for propeller formation depends on the surface density and is independent of 
the optical depth. 

First, we show formation of a propeller in the low surface density model. Figured] shows 
a snapshot of the low surface density model at t = l.OTk, where Tk is the Keplerian period 
27r/0. A propeller-shaped feature is clearly visible in the weak gravitational wakes. The 
surface density decreases considerably in the two lobes downstream of the moonlet. They 
are aligned in the orbital direction and are symmetric about the moonlet. The minimum 
densi ty in the propeller is located at about x = ±2rn = 456 m and y = =f4 km (ISeiB et al. 
20051 ). The radial separation between the two lobes is about 4t"h and the radial width 



of a single lobe is about 2r H rtSpahn fc Sremcevic 2000 ; Sremcevic et al. 2002 ; SeiB et al. 



20051 ; ITiscareno et al.ll2008l ; iLewis & Stewart! 120091 ). The length of the lobe in the azimuthal 
direction is about 5 km. We define the propeller region as the region of —3 < x/rn < — 2 



and 1 5 < y/r-R < 20, which is the typical region of a propeller lo be ( jSpahn fc Sremcevic 
2000l : ISremcevic et~al1l2002[ ISeifi et al II2OO5I ILewis fc Stewart! l2009n . The ratio of the time- 



averaged surface density in the propeller region to the initial surface density E/S is 0.17. 

Next, we demonstrate non-formation of a propeller in the high surface density model. 
Numerical simulations rev eal that a strong wake structur e due to gravitational instability 
forms for this parameter ( !Salolll995t iDaisaka fc Idalll999l ). As expected, until t = 0.5Tk, 
gravitational instability occurs and wakes start forming. Figure |5] shows a snapshot of the 
high surface density model at t = 4.0Tk. Strong gravitational wakes form and no propeller 
structure is clearly observed. The ratio of the time-averaged surface density in the propeller 
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region to the initial surface density S/So is 0.41. The gravitational wakes seem to be almost 
unaffected by the moonlet. The typical distance between wakes is approximated by the 
critical wavel ength of gravitational instability X r r = 47r 2 G£/f2 2 , where £ is the ring surface 
density (e.g., lJulian &: Toomrd Il966l ; ISalol Il995h . In the high surface density model, the 
typical distance between wakes is A cr = 459 m, which is larger than the Hill radius of the 
moonlet. 



3.2. Condition for Propeller Formation 

The numerical simulation results indicate that propeller formation depends on the sur- 
face density. Below, we derive the formation condition for propellers and confirm its validity 
by performing iV-body simulations. 

The clumps in gravitational wakes typically have a mass of ~ £A 2 r . If the clump 
mass is greater than the moonlet mass M, the gravitational wakes may not be affected by 
gravitational scattering due to the moonlet. Comparing the clump mass with the moonlet 
mass, we obtain the following condition for propeller formation 



(2) 



Note that this condition can also be derived by comparing the typi cal velocity due to scat - 
tering by the moonlet with the velocity dispersion of ring particles (ILewis &: Stewartll2009l ). 
If the velocity dispersion of particles exceeds the Hill velocity r^Q, which is a typical velocity 
due to scattering by the moonlet, the gravitational wake should be almost unaffe cted by the 
moonlet. The veloc ity dispersion of ring particles c is determined by Q = 2 (jSalol 11995 
Daisaka fc Idalll999l ). Thus, the ratio of c to r^fl is about A cr /rn. 




Fig. 1. — Snapshot of low surface density model at t = 1.0T K . The initial surface density is 
So = 62g/cm 2 . The circle at the center of the computational box is the moonlet. A propeller 
structure is visible around the moonlet. 
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Using the moonlet radius and the ring surface density, Equation ([2]) can be rewritten as 
\lUn 5 a 6 J 

= 167g/c m * (°-) ( ' Y (^-^Y 3 (-*-) (3) 
' \1.5 J \1.3x 10 s km/ V0-9g/cm 3 / \l00m) KJ 

where C is the non-dimensional constant of order unity (see discussion below) and p and R 
are respectively the density and radius of the moonlet. 

To check the validity of this condition, we perform iV-body simulations for various S 
and R. To investigate a wide range of E, we adopt equal-sized (10 m) ring particles, which 
reduces the number of particles. The density of particles is 0.5g/cm 3 . We assume that the 
restitution coefficient is constant at e = 0.1. The moonlet density is 0.7g/cm 3 . 

We perform 25 simulations with different ring surface densities and moonlet radii. We 
vary the moonlet radius from 60 m to 250 m and the surface density from £ = 133g/cm 2 to 
X = 670g/cm 2 . The ratio of the radius to the Hill radius is 0.64. The simulation time is 
6Tk- We consider that a propeller forms if S/So < 0.2 in the propeller region. We confirmed 
that a propeller-shaped structure is clearly observed when this condition is satisfied. If this 
condition is not satisfied, no distinct steady propeller structure is observed (although the 
ring particles are affected by the moonlet to some extent). Figure [3] summarizes the results. 
It clearly shows that Equation is consistent with the simulation results. We find that 
C = 1.68 explains the S/S = 0.2 boundary well. We also find that C is approximately 
proportional to the boundary value of S/S . 

For a moonlet with a radius of 150 m and a density of 0.9g/cm 3 , the critical surface 
density is 308 g/ cm 2 . In the low surface density model, the surface density is 62 g/ cm 2 , which 
is smaller than the critical value. In this case, a propeller is clearly observed. On the other 
hand, in the high surface density model, since the surface density is 414g/cm 2 , which is 




y ( km ) 

Fig. 2.— Same as for Figure (TJ but for the high surface density model with £ = 414 
g/cm 2 at t — 4.0Tk. Strong gravitational wakes form and a propeller structure is not clearly 
observed. 
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Fig. 3. — Condition for propeller formation in Saturn's B ring in the R-Y, plane. Filled 
squares denote models in which clear propellers form and crosses denote models in which 
no propellers form. We show the ratio of the time-averaged surface density in the propeller 
region to the initial surface density S/E at each point. The solid line indicates the condition 
for propeller formation estimated by Equation ([3]) with C = 1.68. 
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greater than the critical value, a gravitational wake prevails over a propeller. Consequently, 
no propeller is observed. 

Figure H] shows the critical surface density plotted as a function of distance from Saturn 
for moonlets with p = 0.9g/cm 3 and R = 10, 30, 100, and 300 m. Since A cr increases with a 
faster than rn, S cr d ecreases with increasing a. The typical surf ace density in the A ring is 



about 40g/cm 2 (e.g., Esposito et al.lll983t iTiscareno et al.l 120081 1 . Assuming p = 0.9g/cm 



we obtain the critical moonlet radius for propeller formation to be 20 m. Known moonlet s 
with propellers have radii ranging from 20 m to a kilometer ( ITiscareno et al.ll2008l . 120101 ). 
Equation (j3J) is consistent with this observational data. 

Note that this condition is not applicable to rings with low surface densities in which 
gravitational wakes do not develop well. In this case, g ravitational scatt e ring o f individual 
particles is important and the condition discussed in iLewis fc Stewartl ( 120091 ) should be 
applicable, namely that propeller formation is controlled by the ratio of the maximum ring 
particle mass to the moonlet mass. 



4. Summary and Discussion 



We have performed local iV-body simulations to investigate the formation of a propeller 
by a moonlet. By performing real-scale simulations, we demonstrated that in the B ring a 
moonlet with radius R = 150 m forms a propeller in a low surface density ring (60g/cm 2 ), 
whereas it does not in a high surface density ring (414g/cm 2 ) in which the ring dynamics are 
dominated by gravitational wakes. These results indicate that propeller formation depends 
on the surface density. Comparing the moonlet mass with the typical mass of a gravitational 
wake, we derived the condition for propeller formation that the characteristic length of the 
gravitational wakes given by the critical wavelength for gravitational instability be shorter 
than the Hill radius of the moonlet, A cr < rn- We confirmed this by iV-body simulations. In 
a ring with gravitational wakes, the characteristic length of wakes is more important than 
the size of individual ring particles since a wake is a coherent structure. Our condition is 
consistent with observations of propel lers in Saturn's A ring which revealed observational 
signatures of gravitational wakes (e.g., French et al. 20071 ). 



The Ca ssini spacecr aft recently discovered a new putative "moonlet" in the B ring 
(S/2009 SI) JPorcoll2009h . The diameter of S/2009 SI is inferred to be 300 m if it is orbiting 
in the same plane as the ring particles. Its radial distance from the center of Saturn is 



116,914 km, which is at the outer edge of the B ring (jSpitale &: Pored |2010| ). Surprisingly, 
no propeller was observed around S/2009 SI. For S/2009 SI, the critical surface density 
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Fig. 4. — Critical surface density as a function of distance from Saturn for moonlets with 
R = 10 m (solid curve), 30 m (dashed curve), 100 m (short dashed curve), and 300 m (dotted 
curve). The moonlet density is 0.9g/cm 3 . 
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for propeller formation is estimated to be S cr = 308g/cm 2 , which is higher than the typical 
surface density of the outer edge of the B ring. This suggests two possibilities: the surface 
density around S/2009 SI is locally high enough to prevent propeller formation or S/2009 
SI is not a moonlet but a transient feature such as meteoroid impact or temporary clumps. 
Further observation is necessary to determine the nature of S/2009 SI. 

For simplicity, we assumed that the moonlet is fixed. If the moonlet can move, the 
moonlet will b e scattered by gravita tional wakes, which leads to stochastic change in the 



moonlet orbit (ILewis fc Stewartll2009i ). This random motion of the moonlet may hinder pro- 
peller formation. It is important to determine the shape of propeller structures to interpret 
observations. The shape may depend on various ring parameters such as the ring viscosity 
and the size distribution of ring particles. We intend to investigate these problems in a 
subsequent study. 

Observations of propelle r structures reveal bri ght regions rather than low-density fea- 



tures such as a propeller gap (ISremcevic et al.ll2007| ). However, the implied density enhance- 



ment near the moonlet has not been observed in numerical simulations (ILewis fc Stewart 



20091 ). One possible explanation for the br ightness enhancement is collisional debris released 



from ring particles ( ISremcevic et al.l 120071 ). Gravitational scattering of ring particles by the 
moonlet increases the impact velocity between ring particles. Thus, debris may be released 
from ring particles in the vicinity of the moonlet in high-velocity collisions, which increases 
the reflectivity. As shown in Section 13. 2\ the condition for propeller gap formation is equiv- 
alent to that the Hill velocity of the moonlet r^Q is higher than the velocity dispersion of 
ring particles. Therefore, it may be applicable to the condition for local debris enhancement. 
We intend to study this in the future. 

We wish to thank Joseph A. Burns for stimulating discussions on the observations of 
propellers by the Cassini spacecraft. E. K. is grateful to the Isaac Newton Institute for 
Mathematical Sciences for its hospitality during his visit at the beginning of this study. The 
numerical calculations were performed on the GRAPE system at the Center for Computa- 
tional Astrophysics of the National Astronomical Observatory of Japan. 
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